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Abstract 

We evaluate the limit distribution of the maximal excursion of a random walk in 
any dimension for homogeneous environments and for self-similar supports under the 
assumption of spherical symmetry. This distribution is obtained in closed form and is 
an approximation of the exact distribution comparable to that obtained by real space 
renormalization methods. Then we focus on the early time behaviour of this quantity. 
The instantaneous diffusion exponent v n exhibits a systematic overshooting of the long 
time exponent. Exact results are obtained in one dimension up to third order in n -1 / 2 . 
In two dimensions, on a regular lattice and on the Sierpihski gasket we find numerically 
that the analytic scaling v n ~ v + An~ u holds. 

Keywords: random walk, maximal excursion, finite size scaling, enumeration technique, 
Sierpihski gasket 



The random walk (RW) on a lattice has long been studied due to its widespread applications 
in mathematics, physics, chemistry and other research areas. It turns out that despite the huge 
amount of accomplished work, it still remains a thriving research topic. Lots of results can be 
obtained in the continuum limit (Brownian motion) but results for RW on a lattice often yield 
drastically different behaviour - as it is the case for the winding angle distribution |l| - or, at 
least, unusual finite time convergence properties. In the present work we investigate a central 
quantity for RW, the maximal excursion from the origin at time n, M n = max(||x m ||,0 < 
m < n). This random variable is of great interest in many practical purposes such as the 
control of pollution spread, propagation ranges of epidemics, tracer displacement in fluids, 
the radius of gyration of polymer chains |2|, ||, or of lattice animals [|| ||| or other extreme 
statistics. A great deal of work was devoted to the first passage time (FPT) statistics which 
is a closely related quantity. Nevertheless methods used to find the exact FPT distribution 
in one dimensional inhomogeneous environments || [7], |8], [| do not help to get a closed form 
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of the exact distribution of M n . Except for the one dimensional case, only the leading order 
asymptotic expressions (as n — > oo) are available. It was proved long ago by Erdos and Kac 



lOfl , that in this limit the distribution of M n coincides with that of the Brownian motion. This 
result appears as some kind of a central limit theorem. However, it does not deal with centered, 
reduced variables. Moreover it offers no practical access (for physically motivated purposes) to 
the convergence speed towards the limit law. The only global estimates available for M n are the 
laws of iterated logarithm of Khinchine and Chung |TTJ for the one dimensional random walk, 
claiming that although all the distributions have the same limit form, the intrinsic uncertainty 
on M n increases with n. Hence, it is not clear what the finite time behaviour of the maximal 
excursion M n is. It is our aim to clarify this point. 

In this article, we first derive the leading order expression for the distribution of M n in a 
generalized form. This expression is shown to apply also on self-similar structures. Then we 
proceed to the next leading order expansion for short time. In this regime, the first moment 
of M n scales as (M n ) ~ n Un where v n is the effective instantaneous diffusion exponent, and 
tends to u n — > v as n — > oo {y = 1/2 on regular lattices, v = ln2/ln(<i + 3) on the d- 
dimensional Sierpihski gasket). We show numerical evidence that the effective instantaneous 
diffusion exponent v n approaches the limiting value v according to v n — v ~ rT v . This result 
is valid for both regular and self-similar lattices. We finally discuss this point in the context of 
other problems of statistical physics. 

First let us briefly recall the formulae for the maximal excursion of a <i-dimensional Brownian 
motion r t , that is M t = max(||r u || 2 , u < t), where ||r|| 2 = y rf is the Euclidean distance. 
The limit distribution is denoted by P d (a,t) = Pi{M t < a}. The calculation goes through 
the solution of the diffusion equation in d dimensions with spherical symmetry and absorbing 
boundaries on the hypersphere of radius a. Let U(r, t) be the probability density function for 
the position vector r of the walker relative to the origin at time t, without ever crossing the 
hypersphere boundary at distance a. Then U(r, t) satisfies the diffusion equation 

d t U(r,t) = ±V 2 T U(r,t) (1) 

where is the ci-dimensional Laplace operator. The diffusion constant is set as 1 / (2c?) so that 
the solution corresponds to a simple random walk on Z d with a time r between steps and a 
lattice spacing y/r in the limit r —>■ 0. The boundary condition is that U(r, t) = for | |r| | 2 = a, 
and the initial condition is 

U(r ,0) = *(r) = ^W , 

where A d is the surface area of the unit hypersphere in d dimensions and 5 + is the (one-sided) 
delta function. The probability of remaining inside the hypersphere up to time t, P^a, t), is 
the volume integral of U(r, t) over the hypersphere. Due to spherical symmetry, U(r, t) is a 
function of r = ||r|| 2 only, which we now denote by U(r,t), so that, from ([]]) 

d t U(r,t) = -JL^drr^drUfat), (2) 

with U(r,0) = U(a,t) = and 

P d (M)= / A d r d - x U{r,t)dr 
Jo 



2 



The solution of (0) is given in the form of an infinite eigenfunction expansion. This calcu- 
lation can be done for self-similar lattices in the framework of the O'Shaughnessy-Procaccia 
approximation [12|. It consists in assuming a spherical symmetry of a fractal object, and in 
introducing an effective diffusion coefficient D = D r~ 2+1 / U computed from the solution of the 
stationary diffusion problem on self-similar lattices without angular dependence. Thanks to 
this approximation, an analytic approach can be pursued. The final distribution, denoted P^u 
for general v is obtained in the Laplace domain in closed form 



P d .u(a,s) = - 
s 



•>l—dv 



1 - 



dv-l 
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T(du] 



(3) 



Here I n (x) is the modified Bessel function of order n. From this formula all moments are plainly 
computed 



Putting Dq = l/2d and v = 
lattices in one fL0| , two |L3| 
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we easily recover the known distributions on regular 
01 dimensions. A similar method was used in [|14l to 



solve the first passage time problem in the presence of a steady potential flow. It is worthwhile 
mentioning that the limit d — ► oo for v — 1/2 in @ yields Poo(o, s) = - (1 — exp(— a 2 s)) to 
leading order. Hence, as intuitively expected, the maximal excursion of a random walker is 
exactly known in infinite dimension and peaks at a = \fi. 

On self-similar lattices [y ^ 1/2), Equation (|3|) is only an estimate of the exact limit law. We 
have compared distribution (|3|) with the distribution obtained by real space renormalization 
group (RSRG) techniques for the 2D Sierpihski gasket 
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\ [L7j]. Both laws turn out to 
approximate the exact distribution to the same order (see Figure |^ and the discussion in fll8"l ). 
We have also evaluated the moments. To first order, they behave as (M k ) ~ n ku , so we define 
the normalized moments (M h ) = (M k )/n ku which tend to a constant asymptotically. The first 
two normalized moments obtained from (Q), (M) = 1.20 and (M 2 ) = 1.59, should be compared 
with the moments obtained from the RSRG method ((M) = 1.19, (M 2 ) = 1.57) and with the 
numerical estimates ((M) = 1.28, (M 2 ) = 1.84). Both theoretical formulae underestimate the 
actual values |fT§||. This can be understood as follows. Strictly speaking, no limit distribution 



can be defined for M n but, for consecutive time series n, 5n, 5 n, . . ., P n (a/n v ) is left invariant, 
because if a random walker takes T steps to leave a triangle of size R, it needs a time 5T to 
leave a triangle twice bigger. Thus P n (a) fulfills the scaling relation P n (a) = P 5n (2a). However, 
between these times and for fixed a/n u , the rescaled distribution P n (a/n v ) has a log-periodic 
variation. This log-periodic behavior has been known for some time for lattice random walks 
T9| . Both function ([3]) and the RSRG result give the probability to stay in the triangle of size 
R = 2 l , that is one minus the probability to reach sites at distance R + 1 from the origin. As 
observed in Figure (Q) this corresponds to an extremum in the oscillation of P n (a/n v ) (leftmost 
dashed line of Fig. H) rather than to an average. 

Now we turn our attention to the convergence speed towards the asymptotic law (^]). For 
convenience we study the case of the discrete time random walk on the lattice, where analytical 
results can be obtained in one dimension and an exact numerical approach is possible in higher 
dimensions. We focus on the instantaneous diffusion exponent v n which furnishes an information 
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Figure 1: Rescaled distribution of the maximum excursion P n (a/n u ) versus the reduced variable 
r/n v for n = 7500 (x) on the two dimensional Sierpihski gasket, compared to the analytical 
prediction (|), bold line, and to the RSRG result, thin line. Dotted lines: P n {a/n v ) for fixed a 
and varying n: the rightmost curve is computed for a = 2 5 while the leftmost curve corresponds 
to a = 2 5 + 1. Inset A: difference between and the RSRG prediction. Inset B: P n (a/n u ) for 
n = 1000 (A) and n = 1500 (+). The curve for n = 7500 = 5 x 1500 (x) is exactly superposed 
to the curve for n = 1500. 



about the convergence speed of the moments. The numerical estimation of v using a Monte- 
Carlo sampling can lead to false conclusions as in [ED] (see [0]). Here we use exact enumeration 
methods and therefore we avoid such problems. 

The exact solution of the problem in one dimension is obtained by solving the master equa- 
tion with absorbing boundaries at points ±a with the use of a Fourier development (obtained 
22 with a minor misprint), but the moments cannot be calculated in a straightforward 
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manner from this expression beyond first order. We derived another form of the distribution 
by a recursive use of the reflection theorem. The probability density of the maximal excursion 
at step n reads 

oo 2k 

P n {M n = a) = 2 ^(-l) fc [ Pn ({2k + l)a) + p n ({2k + l)(a + 1)) + 2 ^ Pn ({2k + l)a + i)] (5) 

k=0 i=l 

where p n {x) is the probability density function for the discrete random walk to be at position x 
(which is null for |x| > n). Formula (^) allows a convergent expansion of the first moments in 
powers of n -1 / 2 , which exist, since P n {M n = a) is an analytic function of n l l 2 . At order n" 1 ^ 2 , 
divergent series are encountered which can be summed by classical methods [23||, yielding 



/,tx F™ i i 2n _/r. 
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<M»=2a.-jf + £±i + 0(-jg| (7) 



where G = 0.9166 ... is the Catalan constant. In the calculation of the second cumulant the 
terms of order y/n and 1/y/n cancel, as expected. The distribution Pt{M t = a) for continuous 
time random walks (CTRW) follows plainly from since Pf(M t = a) = P n (M n = a)U t (n), 
where n 4 (n) is the probability for the CTRW to perform n steps in time t. Numerically, a 
series expansion similar to (0) is found for an exponential distribution of waiting times. The 
exponential distribution is particular because it is the only one for which a master equation 
formulation and a CTRW on the same lattice are isomorphic ||24| . A striking feature of the 



random variable M compared to other extremal quantities at finite times is that the leading 
order expansion of (M*) scales as n k l 2 + cte.n^ -1 ^ 2 and not as n k ^ 2 + cte.n fc//2_1 , hence finite 
size effects persist for a large number of steps. 

In two dimensions, no exact result is available for finite times and the analyticity of the 
probability density P n (M n = r) is not obvious. Hence we investigate this case numerically. It is 
possible to perform an exact enumeration of the walks by studying the joint probability density 
of the position and maximal excursion P n (x,y, M) on the square lattice Z 2 . We can compute 
P n (x, y, M) in the region < x < M, < y < x only, due to symmetries. We use the family of 
metrics d p (x) = (^2 i (\xi\ p )) 1 ^ P to compute the maximal excursion from the origin of the lattice. 
In Figure ^| we plot the instantaneous exponent v n with three classical choices of metric: d\, 
g?2 (Euclidean distance) and d^ (max distance). The metric d\ and both induce a strong 
overshooting of v n with the limit value 1/2. The curves have the same shape as that of the 
one dimensional case, also potted for reference in Figure |2|. The first two moments have many 
features in common with their ID equivalents. For example, using the metric d^ we find that 
the series expansion (M*) = J2^=o m k- P (Vn) k ~ p holds for both first and second moments up 
to third order, with mj = -0.50, = 0.322 and m\ = -1.083, m 2 = 0.95, m\ = -0.33. 
The leading order terms are exactly evaluated from the asymptotic results and read m\ = 
1.0830, m\ = 1.3048. In the Euclidean metric d 2 , we find a drastic change in the shape of the 
curve v n (Figure |2|). The instantaneous exponent approaches \ from below and remains below 
| at time n = 400. This phenomenon is a lattice effect. We show this fact by computing u n in 
an off-lattice random walk model (Figure 0). Since it is not possible to use exact enumeration 
techniques in this case, we resort to a Monte-Carlo simulation. We inspect 2.10 s random walks 



with fixed distance increments and a uniform distribution of the angles (Pearson walks p2|). 
In this situation, v n in metric d 2 is very close to that obtained in metric d\ and c^o, and it 
decreases towards |. Both lattice and off- lattice models should give equivalent results once the 
discretization effects are smoothed out. Therefore v n should ultimately approach | from above 
in the on-lattice model. We have investigated the change of v n when varying continuously the 
metric d p with 1 < p < oo on the lattice (Figure 0). For large enough values of p (p > 50), we 
do observe that the curve crosses the value ^. In the metric d 2 , however, the time needed for 
v n to cross | should be enormous. 

This result shows that the definition of the metric strongly influences the convergence prop- 
erties of the maximal excursion on regular lattices. The two natural metrics for the square 
lattice, d\ and d^, lead to a behaviour similar to that observed in the continuum model. 

On the Sierpihski gasket we enumerate the walks starting from the top of the biggest 
triangle up to a fixed time. The metric chosen here is the chemical distance from the origin. 
For each maximal excursion r we compute the probability of remaining below r after n steps, 
P^{r). Unlike the transfer matrix method, this method works only at fixed time, but allows 
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Figure 2: Figure 2: Instantaneous exponent v n (averaged over two consecutive steps) versus 
number of steps on the square lattice at two dimensions. Enumeration up to step n = 400 in 
metric d\ (bold solid line) and metric d^ (solid line). The one dimensional situation (caliper 
diameter) is given for reference (dot-dashed line). v n is also computed using the metric c?2 
(Euclidean distance, dashed line), di.05 and d 5Q (solid lines). 



us to discard the long transient regime and to compare the exact limit distribution to its 
spherically symmetric approximation. We have computed the instantaneous exponent v n up to 
step n = 10 4 on a gasket of size 256 (cf. Figure [|). Like the moments, v n displays a log-periodic 
oscillation persisting in the long time regime with an amplitude less than 8 10~ 3 . v n tends to 
the asymptotic value v — for long time. It seems that on a very general class of lattices 
the finite time behaviour of v n and therefore of (M^) is an analytic function of rT v . This fact 
lacks a clear physical understanding. The real space renormalization results do show that n u is 
the well defined time scale for this problem but the exact evaluation of finite size effects is not 
accessible from this method. To assess this hypothesis we have smoothed out the log-periodic 
oscillations of v n . For a log-periodic function f(x) = f(Tx), one can define z = ln(x) and 
f(z) = f(x) so that the running logarithmic average reads 

Jz Jx V 

Using a discrete form of this average we write 

5n 

1 P i 

V5n = = > 

2—ii=n i i=n 

We tried to fit V n using 

V n = v + An~ v + Bn~ 2u + o(n" 2v ) (8) 

and we found a very good regression for A = 0.082±0.002 and B = 0.18±0.1. However, the best 
fit with only one power law is u + cte.n~ a with a = 0.49, so that strictly speaking a nonanalytic 
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Figure 3: Instantaneous exponent v n versus number of steps in the off-lattice model at two 
dimensions. Monte-Carlo simulation up to step n = 100 in metric d% (bold solid line) d<i (dashed 
line), and doo (solid line). The results for metric d\ and d^ are almost indistinguishable, as 
expected. Inset: Log-log plot of V n — v which shows the n" 1 / 2 scaling. 



short time dependence cannot be ruled out. The underlying assumption in the computation of 
V n is that the regular log-oscillatory pattern of v n is additive. This assumption does not hold 
because the averaged exponent V n still shows some oscillation. The local exponent a fluctuates 
between 0.40 and 0.52, which does not allow to confirm unambiguously the hypothesis of the 
analytic behavior of function of n u . 

We would like to point out that in the context of lattice animals, the quantity (M%), or 
equivalently the "caliper diameter" (average spanning diameter of lattice animals once projected 
on a fixed axis) displays a sub- leading order behaviour which scales as n^ k ^ v [|j], aside from the 
well-known non analytic subleading term, and can be interpreted as a 'surface contribution'. 
In the case of the maximal excursion of a random walk, we have proved in one dimension 
and evidenced through enumerations in higher dimensions that the early time instantaneous 
exponent v n is systematically above its limit value v with a leading order development v n ~ 
v + Arr v where A depends on the precise choice of the metric. This result is consistent with the 
fact that the corrective scaling to the moments due to finite size effects includes only terms of 
the form (n~ u ) p , peN, which was proved in one dimension and which can also be interpreted 
as a surface contribution. 

In conclusion, besides its intrinsic interest, the maximal excursion of a random walk shares 
difficulties which are often encountered in physical problems dealing with finite size series. In 
certain cases, power law exponents inferred from the finite size series expansions should be 
considered with caution, as might be the case for directed percolation series. 
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Figure 4: Convergence of the instantaneous exponent for the first moment v n (bold line) towards 
the limit value In 2/ In 5 (dashed line) on the two dimensional Sierpihski gasket. The running 
average V n is also plotted (dashed bold line). Inset: Log-log plot of V n — v and fit 
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